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Magnetic fields appear to be a generic feature of the early universe and are a natural source of secondary 

CMB non-Gaussianity. In recent years the statistical nature of the stresses of a primordial magnetic field has 

been well studied. In this paper we confirm and extend these studies at one- and two-point level, and present 

analytical results for a wide range of power-law spectra. We also consider two non-power law cases of interest: 

a blue spectrum with an extended damping tail on small scales, which could be generated by the non-linear 

mixing of density and vorticity; and a red spectrum with a damping tail on large scales. We then briefly consider 

f^ the CMB impacts that result from such fields. While this paper focuses on the one- and two-point moments, 

y—{ the techniques we employ are designed to ease the analysis of the full bispectra induced by primordial magnetic 

O fields. 
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I. INTRODUCTION 



Large-scale magnetic fields appear to be inevitable. Fields have been observed on galactic and cluster scales and likely 
exist on supercluster scales 1 1 -4|. There are a variety of potential generation mechanisms in both the early and late universe, 
many of which may well have been in operation, ranging from so-called Biermann batteries (e.g. 15] ID) to fields generated at 
reionisation It HToI . The number of well-motivated magnetogenesis mechanisms would appear to make the existence of magnetic 
\^ fields across a large range of scales inevitable. In particular, plasma processes in the early, pre-recombination universe 1 1 1 - 17| 
wherein linear scalar perturbations drive nonlinear vorticity, producing currents in the proton/electron plasmas, would suggest 
O , that at least some magnetic fields before recombination are inevitable even if the detailed predictions in the various models dififer 
I somewhat. The fields extent in the present universe need not necessarily be primordial; the fields observed on cluster scales can 

be accounted for in principle by non-linear processes before recombination, by production at reionisation or at recombination 
||T81[T9J , and it seems probable that each of these processes have contributed. There also many exist potential mechanisms in the 
jrt extremely early universe, such as magnetogenesis at a phase transition (for example Il20l424l '). 

I— I However, all of these fields will have a distinctive, blue power spectrum t25), and the most recent bounds on fields present 

at or before recombination rather favour a red spectrum |26, 27|. It seems reasonable, therefore, to also consider mechanisms 

►^ that could produce such red fields. To achieve a near-scale invariant spectral index these will necessarily be produced during an 

p^ inflationary epoch [251, and there is a large variety of potential genesis mechanisms that generate a field from the breaking of 

QQ the conformal invariance of the electromagnetic field, with |28-33| being a few examples. 

0^ Constraints on a field produced during the inflationary epoch have naturally focused on the cosmic microwave background 

C^ (CMB) ( Il26ll27ll34l - l36l . for example) and have typically found that the field strength today is at most of the order of nanogauss. 

ly-N A stochastic background of primordial magnetic fields, such as would be reasonably expected from an inflationary process, 

(^ impacts on the CMB in two chief ways. Firstly the magnetic field itself acts as an additional radiative fluid which, should there 

(^^ be enough energy in the field, impacts on the expansion rate of the universe. However, by far the most significant impact comes 

^~~^ from the interactions of the magnetic field and the ionised plasmas extent before recombination, and the scattering between the 

^ magnetic fields and the geometry. There have been numerous studies into magnetised cosmological perturbation theory; for 

some interesting studies see for example |3 , 37-47 1 and their references. Clearly the presence of the magnetic field will impact 

on the CMB angular power spectra. While the impact on the temperature auto-correlation and on the temperature/polarisation 

t^ cross-correlation will be small except at larger multipoles, a magnetic field also excites B mode polarisation. The B mode angular 

power spectrum would therefore be of great interest for the study of primordial magnetic fields. See for example |40, 48-581 

and their references for studies of the magnetised CMB angular power spectra. 

The B mode spectrum produced by a primordial magnetic field is, however, likely to be nearly degenerate with the gravitational 
wave spectrum produced by inflationary models ll52l and even with forthcoming results from for example the Planck and QUIET 
experiments it is possible that we will be unable to disentangle the impact of magnetic fields from other contributions and from 
foreground signals. An alternative is provided by the non-Gaussianity of the magnetic field, studied in for example 145. .59.-61 1 ; 
since the stress-energy tensor and the Lorentz forces generated by a magnetic field are non-linear, the statistics imparted on the 
perturbations by the magnetic field are necessarily non-Gaussian, regardless of the statistics of the underlying magnetic field. 
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Furthermore, we would expect the greatest impact to be on large scales, observations of which with the WMAP satellite are 
already cosmic -variance limited |62|. Study of the magnetised CMB bispectrum is, however, unfortunately a difficult issue. 
It requires knowledge of the full 3-point moment of the stress tensor of the magnetic field (the "shape function" or "intrinsic 
bispectrum") and to date only a few particular configurations have been studied. While these cases are probably sufficient for 
nearly scale-invariant magnetic fields, for which power appears to be concentrated on the so-called local configuration (with the 
bispectrum geometry set by three wavevectors ki » -k2 and k^ ^ 0), it has not been possible to fully test this. Furthermore, for 
the blue cases expected to arise in non-inflationary magnetogenesis models neither the local, nor indeed any other, configuration 
will necessarily dominate and knowledge of the full intrinsic bispectrum will be required. 

To consider the impact of a magnetic field on the CMB, then, we require a clear understanding of the statistical nature of 
the stress tensor (and, equivalently, the Lorentz forces). In fSOl and p^"] (hereafter BC05 and B06) we studied the 1-, 2- and a 
special case of the 3-point moments of the stress tensor of a Gaussian magnetic field with a power law power spectrum with index 
hb, employing a mixture of realisations of the magnetic fields and numerical integration. The nature of the 2-point moments 
changes depending on whether hb > -3/2 or hb < -3/2 and we considered two cases in particular For a field with a flat 
spectrum we found the 2-point moments with both methods, demonstrating good agreement and describing the gross features 
of such "ultra-violet" fields. For a field with a strongly red spectrum we only employed realisations and could only draw the 
broadest of conclusions. In B06 we further considered a "damped causal" power spectrum, resembling that produced by plasma 
processes in the pre-recombination universe, and the statistics of the Lorentz forces (as subsequently studied in, for example, 
iSTi). More recent studies have improved the approach to numerical integration |55 1, and extended it to a wider range of spectral 
indices; and have even found analytical solutions for particular cases |56 57 1. We argue in this paper that the most important 
spectral indices are at ng = 2, hb - 0, n^ = -3/2 and hb - -5/2, all of which are included in their solutions. 

However, while this means that we now have solutions from analysis, numerical integration and from statistical realisations 
work remains to be done. In principle the realisations are the most general approach, and can be applied to any field, possessing 
any statistical nature and any power spectrum. However, they are severely compromised by the mode coverage on large scales 
associated with the finite size of the grid. The extent of this issue was not considered before. Due to its reliance on Wick's 
theorem, numerical integration can at present only be applied to fields with Gaussian statistics, although arbitrary power spectra 
can be employed. Our approach in BC05 and B06 was restricted to blue power-law spectra with n^ > -3/2 which is overly 
restrictive. The approach of [5 5 J is more general but relatively unwieldy and has been applied only to power-law spectra. Finally, 
while 1571 found analytic solutions for spectral indices of interest, they are missing in particular the case hb - -2. As there 
are only two analytic solutions in the region riB e (-3, -3/2), at hb - -5/2 and hb - -2, this is an important omission. The 
subsequent approach to CMB observables, given these spectra, has been rather piecemeal and model-dependant. 

It is our intention in this paper to reconsider many of these issues, always focussing on retaining as much generality and 
flexibility as possible. In particular, we wish to explicitly decouple the consideration of the the intrinsic magnetic statistics from 
the evaluation of the CMB angular power spectra; in this manner, one may build the statistical nature of the magnetic stress 
tensor in any manner one wishes, and wrap it onto the CMB employing independently-derived analytical or numerical transfer 
functions. Further, we wish to approach the issue of the intrinsic statistics in an extensible manner, unifying to some degree the 
approach to analysis and to numerical integration. In doing so we extend the set of useful analytic solutions, finding solutions for 
all hb between «b = 3 and ns = -5/2. While we will not generalise beyond underlying Gaussian statistics for the magnetic field, 
we also consider two non-power law power spectra of interest. The first is a damped causal spectrum similar to that considered 
in B06, and the second resembles a strongly red field produced during an inflationary epoch with an extended decaying tail on 
large scales. We find analytical solutions for the former and rely on numerical integration for the latter. We then consider the 
CMB, employing a clearly-defined approximation based on the generic form of a magnetised transfer function. In this paper 
we consider for brevity only the 1- and 2-point moments but we emphasise throughout the paper that we have constructed our 
techniques to be immediately applicable to the study of magnetic bispectra, to which we will turn in a forthcoming study. 

We start in section|ll]by introducing our toy model of primordial magnetic fields and briefly discussing the normalisation of 
their power spectrum. In section [nl] we briefly consider the one-point statistics of the stress tensor before considering in some- 



what more detail the two-point moments in section IV In section V A we present our general approach to the magnetised CMB 
angular power spectra, and in section jVBJ apply this to the CMB temperature auto-correlation generated by tensor perturbations. 
We briefly conclude in sectionjVl] Appendix [A| presents those analytical solutions not previously presented in |57|, appendix 
iBlpresents the forms on the largest scales, and appendix Obriefly addresses a comment in 1.5 5 J concerning our previous work. 
Latin indices refer to spatial dimensions and are raised and lowered with the Minkowski metric rjij = diag(l, 1, 1). 



II. TANGLED MAGNETIC FIELDS 



At the linear level a large-scale primordial magnetic field baik, rf) decays as a^{T]). We work with the scaled field Ba(k) = 
a^{Ti)ba(k,Ti), which is constant over time If37ll38l . and assume the electric fields to be negligible. Magnetic fields are damped 
on small scales, primarily from radiation viscosity, and we approximately model this effect with a time-dependent cut-off scale 
kcirj) above which they have no power If the magnetic field is generated at some epoch //in with some cut-off scale A;^'™ and the 



viscous damping is denoted by koirj), the effective damping scale of the magnetic field will then be 

M77) = min(^P™M'?)). (1) 

The magnetic fields are then static for a brief period before the viscous damping scale evolves through the primordial, and 
the field begins to damp on increasingly large scales. Note that there is then a time-dependence in the magnetic field Bfl(x) 
associated with the damping scale; in this section we will write this explicitly. The damping scale freezes at photon decoupling; 
hence ^,.(770) ~ ^c(?7rec)- 

Take as a toy model scaled magnetic fields - not necessarily Gaussian - with the power spectrum 

<B„(k, k,{ri))Bl{\.' , kAn))} = pB(k)Pab(y^)H(kATi) - k)(2nf6(k - k') (2) 

where 

Pab(y^) ^Sah-^ (3) 

is a projection tensor that projects an object onto a plane orthogonal to the wavevector k, VBik) is the magnetic power spectrum 
and in the interests of simplicity we have neglected an antisymmetric, helical component. H(x) is the Heaviside function, which 
we define to vanish for x < and to be unity for x > 0. For a magnetic field with a power-law spectrum 

rB{k)^ABk"\ (4) 

we can immediately see that analyticity of VBPah requires hb > 2, corresponding to magnetic fields generated by causal pro- 
cesses. Fields produced by the generation of second-order vorticity from linear scalar perturbations, as considered with varying 
results in fTP-ri?], will naturally be of this type on the largest scales, with a decaying tail on smaller scales. In this paper we 
model cases both with and without an extended damping tail. 

It is common in the literature to normalise the magnetic spectrum to the average magnetic energy density at the present epoch 
and on a scale A, typically of the order of a megaparsec and associated with cluster scales. Smoothing the field with a Gaussian 
filter f;^{k) - exp(-/l^A;2), one finds (e.g. l|6Tl) that 

B",-^ I rB{k)e-' ' k^dk. (5) 

^ Jo 

For the power-law spectrum (J4]) this becomes a standard integral and has the general solution 



where r{a,x) is the incomplete Gamma function. (Note that this differs slightly from that presented in |56| due a different 
definition of the smoothing function.) If kdrio) » k^, or if n^ —^ -3, this reduces to 



?2 ^ ^Ji 1 ^lnB + 3 



5l-;^^^rpi— , (7) 



in agreement with previous estimates in, for example, pO). From this we can see that we require hb > -3 to keep the integral 
finite. (See for example |25 58 61 , 63] for some further discussion of these issues.) Since the damping scale is typically rather 
small, we assume (J7]i to hold for the power-law fields we consider In cases where an analytic solution to dSll does not exist the 
equation can in principle be solved numerically or approximately. 

An alternative that has recently been employed l.56-,- .58. ^61J is instead to normalise the magnetic spectrum with the mean 
square field itself. 
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PB{k)k^dk. (8) 



The two normalisations are clearly closely related, and in this paper we employ the former. 

In addition to the minimally causal case n^ - 2, we can identify two further power-law spectra of particular interest. In 
BC05 and B06 we studied in detail the "flat" or white-noise field with ng - 0, which highlights the impact of each mode on the 
statistics. A field with hb - -5/2 is also of special interest as it is a half-integer - which will shortly prove extremely useful - 
that lies in the centre of the current bounds on the magnetic spectral index ll26ll27]l and could be generated during an inflationary 
epoch. A field with hb - -2 also satisfies the current bounds, but constraints are weakest for hb — » -3 (e.g. Il63ll ). We typically 
refer to fields with «b > -3 as "inflationary" fields, although in principle such a field can take any hb > -3. 



The magnetic field contributes to the Euler and Einstein equations through the stress-energy tensor 7^(k). As the Poynting 
vector vanishes to first order and the magnetic energy density is equivalent to the isotropic pressure it is sufficient to consider the 
stress tensor T^(k), 



r°(k) = f;(k)5« - f;;(k), where f,,(k) = f B„(k')Bi,(k - k')d^k' 



(9) 



is the self-convolution of the magnetic field. The stress tensor can be separated into the isotropic pressure (scalar trace), 
anisotropic pressure (traceless scalar), vorticity and transverse-traceless (TT) tensor components, which we denote by t, t^, 
T^ and tJ,^ respectively. For more details see B06. The statistics of the impact of the magnetic field can be characterised by 
studying the rotationally-invariant statistics of these components. We consider here the 1-point moments (the probability dis- 
tribution functions, skewnesses and kurtoses of the scalar pressures) and on the 2-point moments (the auto-correlations and the 
cross-correlation ( rrj )). The magnetic field also contributes through the Lorentz force and this has been studied in B06 and in 
Il27ll57l : however, here we focus on the statistics of the stresses themselves. 

Assuming the fields to be Gaussian in nature allows us to both analytically and numerically integrate the expressions for the 
stress power spectra, as done numerically in BC05, B06 and |55|, and analytically as in |56, 57|, and we extend these approaches 
both to a wide range of spectral indices ng and to two non-power law spectra of interest. Some of our analytical results were 
previously known |57|, but others, including the damped causal spectrum, are original to this paper, as are the numerical results 
for the IR-controlled field. 

To aid the analysis, particularly at the one-point level, we also employ extensions of the codes constructed in BC05 and B06, 
creating static realizations of B„(k) numerically. Details of the construction can be found in BC05 and B06, and the grid is 
characterised by its side-length /dim- In earlier studies we employed /dim - 192, while in this paper we will instead employ 
^dim - 256 or, in some cases /dim = 512. The Nyquist frequency of the grid is A:Nyquist = 'dim/2; due to the quadratic nature of 
the stress tensor, we should then take the damping scale to be kc < /dim/4 to avoid aliasing power from large to small scales. In 
Figure [T] we present a sample Gaussian realisation of a magnetic field with spectral index hb = -5/2 along with its isotropic 
and anisotropic pressures, t and ts respectively. All errors presented on quantities derived from statistical realisations are 1-cr 



III. ONE-POINT MOMENTS 

In BC05 and B06 we considered the probability distribution functions of the isotropic and anisotropic pressures of the 
magnetic field, and their skewnesses and kurtoses. The central moments of a distribution P(x) with mean ju'j are defined by 

/i„ = ((x - yUj)" ). The second central moment is the variance, while the (normalised) skewness and kurtosis are 

ri = ^7^, 72 = -^-3 (10) 

A'2 ^2 

where the kurtosis is defined such that yi = 72 = for a Gaussian distribution. 

The isotropic pressure, t oc BjB' of a Gaussian magnetic field is expected to follow a x^ distribution with three degrees of 
freedom. The skewness and kurtosis of ax^ distribution with /? = 3 degrees of freedom are 

12 
72 = — =4. (11) 

P 

In BC05 and B06 we confirmed the x^ nature of the isotropic pressure employing twenty realisations at /dim - 192, finding 
y[ - 1.63 ± 0.01 and yj - 3-99 + 0.05 for ng = 0. The anisotropic stress is harder to characterise but we argued that the 
distribution function should be composed of a mixture of x^ and modified Bessel distributions and should therefore have a 
relatively small skewness. We found numerically that the statistics are dependent on the spectral index, with y|' = -0.24 + 0.003 
and yI' = 1.10 ± 0.01 for «b = 0, and y]' = 0.38 + 0.01 and yj' = 0.86 + 0.02 for «b = -2.9. 

However, we did not probe this dependence in detail and we briefly consider it here. For ng > -3/2 we employ twenty 
realisations at /dim - 256 while for n^ < -3/2 we use twenty realisations at /dim - 512. In the left panel of Figure |2] we plot 
the distribution functions for t and tj . The isotropic pressure is invariant under changing ng to a good degree and we present it 
only for riB - 0, while we plot tj for both n^ = and ng - -2.9. The x^ nature of the isotropic pressure is very clear, as is the 
change in the skewness of the anisotropic pressure between the two cases. 

In the right panel of Figure|2]we plot for a wide range of «b the skewnesses (black) and kurtoses (green) of the scalar pressures. 
The statistics of the isotropic pressure are invariant across the range of rtg, with y[ « 1.633 and -y^ ~ 4. The errors grow more 
significant and the mean of yj reduces as hb — > -3 due to the lack of modes with appreciable power in the most steeply-tilted 




cases. Explicitly, at ng = 0, rl = 1-634 + 0.004 and y^ = 4.011 + 0.031, while at «b = -2.9, y] = 1.620 + 0.084 and 
yl = 3.932 + 0.284. 

The anisotropic pressure behaves quite differently. The kurtosis is invariant at y^^ ^ 1.1 for ng € [-3/2,2], but after a 
transition at ng = -3/2 the kurtosis climbs towards y^^ < 2 for n^ — > -3. The skewness depends even more strongly on n^, 
rising monotonically for all hb as ng decreases. At «g = -3/2 the skewness passes through zero. The magnetic field with 
hb = -3/2 is therefore highlighted as a field of particular interest; for this field, the statistics of the anisotropic stress change 
their behaviour Explicitly, at hb = 0, y\' = -0.238 + 0.004 and y^' = 1-108 ± 0.01 1, and at «b = -2.9, y]' = 0.697 + 0.067 and 
yl' = 1.619 ±0.105. 

Our findings for the isotropic pressure and the anisotropic pressure at n^ = are in full agreement with those in BC05 and 
B06. However, the statistics of the anisotropic pressure at hb - -2.9 reveal a tension. This is an indication of a strong grid 
dependence as one tends n^ towards scale-invariance. Having repeated the analysis with varying grid-sizes across a range of 
spectral indices we can determine that increasing the grid resolution tends to increase the skewness and kurtosis of the anisotropic 
pressure. While this is still the case for /dim - 512, it is to be expected that the true results are close to, but slightly larger, than 
those presented in this paper. 

Perhaps the most important result from this section then still holds - that the field with n^ - -3/2 is selected as one of special 
interest and that at this point the skewness changes sign. This field also remains within observational bounds Il26ll27l . The other 
point of interest we can note is that, as hb — > -3 as seems natural, the skewness of the anisotropic pressure tends towards y[^ < 1 
and the kurtosis tends towards yj' ^ 2. A firmer statement would require further analytical study. 

IV. INTRINSIC POWER SPECTRA 

In BC05 and B06 we derived expressions for the correlations between components of the magnetic stress tensor, which take 
the form of an integral across two power spectra.' However, due to the presence of poles in the integrations for hb < -3/2 we 
only integrated these equations numerically for hb - and relied on statistical realisations for ub < -Jill. These results were 
then compromised by limits on the dynamic range of the realisations and could neither be extended into the low-A: region nor to 
k> kc- In this section we improve our analysis, rederiving and extending work by 1.55. ,57J . 

In BC05 and B06 we pointed out that at the two-point level magnetic fields can be separated into "ultra-violet" and "infra- 
red" fields. Ultra-violet fields are those with hb > -3/2; the stress spectra are dominated by the damping scale, and are 
white-noise on large scales. Infra-red fields have hb < -3/2 and on large scales the stresses obey a power-law spectrum. 
The field with n^ - -lill is therefore also of interest for two-point moments. In Il56ll57l analytic solutions were found for 
Hb e {3,2, 1,0,-1,-3/2,-5/2). These solutions encompass the fields we have highlighted as of most interest; the minimally 
causal case «b - 2, the white noise field ub - 0, the "transitionary" field hb - -3/2, and an inflationary field with hb - -5/2. 

However, further analytic solutions can be found, and the integrations can be re-expressed in a form more explicitly amenable 
to numerical integration. Our approach to these allows us to consider a wide range of power spectra, not necessarily power-law; 
and, most importantly, the techniques presented in this section are readily adaptable to the study of magnetic bispectra, the 
integrands of which are significantly more complicated than those for the power spectra, and which are difficult to simulate 
statistically (see B06). In this section the amplitude of the power spectrum is assumed to be Ag s 1 unless otherwise stated. 

If Ti and T2 denote two components of the stress tensor (isotropic pressure, anisotropic pressure, vorticity or TT tensor), then 
the equal-time correlation between the two for a Gaussian magnetic field is B31 l55l l57l l59l 1611 

<Ti(k)T;(p)) = rx2{k)S(\i - p) = 5(k - p) ivBik'WB (|k - k'l) Ti2{\i, k')d^k'. (12) 

Three angles determine the relations between the wavevectors k and k', 

7 = k-k', A'^k'-k^', yS = k-k'^' (13) 



' In contradiction to a statement in |55|, tliese expressions were exact and did not neglect any terms. See Appendix C for further details. 



and the angular terms in the integrations are given by 

T,, = ^(l+A^'), (14) 

Tr,r, = \{\ + y^-l,[j'+li'-)-(,ypifi + 9y'f), (15) 

Trr, = ^(//2 + 3(r'+/32)-3r///?-2), (16) 

5>v,v = l-2ry + y;u/3, (17) 

Solving this integral is problematic due to the poles that can appear in fBiK)^ ^B(|k - k'|) and Tn- To control this it is useful 
to express k! , k and p in units of the cut-off scale kc, writing k' = akc to ease the notation, and change variables such that we 
integrate across y = |k - ap. This ensures that the poles in !P(|k - a|) can appear only at the lower limit of integration and are 



easier to control numerically.- Equation ( 12 1 can then be recast as 

nk^ r' rik+aj- 



rn{k)^-^\ aVBia) \ n(>6^)^i2(^,«,y)dydfl. (19) 

!<■ Ja=0 Jy=(k-af 



In terms of these units, the angles defined in ( 13 1 are given by 



k^ + a^ -y ky - a k^ - a^ -y k - ay k^ - a^ + y 

2ak ' -v/y 2ay/y ' yjy 2k ^^ 

For convenience, we define the scaled spectrum 

QmQC) = ■^'Pnik). (21) 

The power spectrum VBik) vanishes for arguments less than zero or greater than unity, which restricts the integration limits to 
particular values. These vary depending on the value of k: 

• For k > 2, the lower limit on the integral across y is always greater than unity, and the integral vanishes. Therefore 
Qnik > 2) = 0, as has previously been found numerically (BOS and BC06, ll55l ). and analytically lISTllFTI and is expected 
from the quadratic nature of the stress tensor 

• For k € [1,2), the lower limit can be less than unity, and so there is a non-vanishing contribution to the power spectrum. 
This contribution decreases as A; — > 2 and the range of integration vanishes. If a > A: - 1 then {k — of- < 1; however, 
{k -H fl)^ > 1 at all times. Therefore 

Qn{ke[\,2))^ { aVBia) { 'PB(^[y)Tx2{k,a,y)dy^a (22) 

Jk-l J(k-a)- 

and this tends smoothly to zero as k tends to twice the damping scale. 

• For k e (0, 1), the integral across a can be separated into the regions a e [0,1 - k] and a e [1 - A;, 1]. In the former region, 
(k - a)^ and {k + of- are both always less than unity, while in the latter (k + a)^ > 1. The integral is therefore 



Jn^-k Mk+dy pi pi 

ana) 'P(y/y)ri2(k,a,y)dyda+ aPia) n^Wn{k,a,y)Ayda. (23) 

J(k-a)- Ji-k Jik-a)- 



-1-i r-(k+a)- /-I /-I 

'P(^)Tx2{k,a,y)Ayda+ aPia) 

l(k-a)- Ji-k J(k-a)- 

In the interests of numerical stability, this integral could be further separated depending on whether a < k or a > k. 
• Finally, for k -Q the domain of the integration across y vanishes, and Qnik = 0) = 0. 



2 



This change of variables is similar to that employed in, for example, 1 64 1 where the focus was on the gravitational waves produced by first-order scalar 
perturbations. Such a change of variables might also be useful in studies of inflationary bispectra where again integrals across two or more power spectra 
must be taken. 



We can solve these integrals analytically for all integer and half-integer ng e [-5/2, 3], adding hb e {5/2, 3/2, 1/2,-1/2, -2) 
to the solutions given in |57|, and verify these results with numerical integrations. Most importantly, this means that we now 
possess analytic solutions for both tig = -2 and ng = -5/2, which are infra-red and close to current observational bounds. 
Full analytic solutions not given in [57 1 are presented in Appendix [A] Numerical integration was performed with Monte-Carlo 
routines, using with 4,000-20,000 samples per integration, depending on ng-^ 

In the left panel of Figure [3] we present the stress power spectra as found for n^ = from the analytic solution and by 
employing statistical realisations. The results from the numeric integration are indistinguishable from the analytical solutions. 
While in BC05 and B06 we employed twenty statistical realisations on a grid of size Zdim = 192, here we present instead the 
results from averaging one hundred realisations generated on a grid with /dim = 256. The plots clearly demonstrate that the stress 
spectra are dominated by the damping scale - at each k both the gradient and the amplitude of the field are strongly dependant 
on it. Since this scale evolves in time, this implies that the statistics imprinted by a cosmological magnetic field vary across time. 
We refer to this as "decoherence" and it implies that studying the imprint of such fields on the CMB may require careful study, 
as for defect mechanisms (see for example [65 j). This behaviour holds for all ng > -3/2. For k < kcoh where A;coh ~ A;c/100, the 
stress spectra tend towards the white noise solutions discussed in the earlier literature (e.g. 1401 l63l ). If the viscous damping is 
on a small enough scale these fields can be treated as white noise and the CMB signals recovered in the usual manner However, 
this should be carefully verified to avoid wrongly estimating the impact on the CMB. In contrast to the auto-correlations, for 
^ < ^Coh the cross correlation !P„ tends to vanish and is only significant atk^ k^. 

In the right panel we present the power spectra that are produced by an inflationary field with ng = -512 - this could be 
compared with Figure 4 of BC05. In BC05 and B06 we employed realisations at /dim = 192. Furthermore, the damping scale 
was set at kc - /dim and we relied on the steeply tilted power spectrum to control the aliasing of power from large to small scales. 
This limited us to the study only of A: <s kc- Both of these issues are addressed here; we use /dim = 256 and set kc - /dim/4 and 
employ the average of 100 realisations to reduce the errors further. Agreement with the analytic solutions for A; > A;^. is extremely 
good. However, for k < kc we remain compromised by an infra-red damping arising from the infra-red cut-ofF associated with 
the finite grid size. 

In contrast to the ultra-violet case, the gradients of the spectra produced by an infra-red field become independent of the 
damping scale at A: < kcoh where kcoh - kc- They also scale as Vnik) oc k^""^^, in agreement with previous papers. The statistics 
imprinted by such a field on relatively large scales therefore do not evolve through time, and the CMB signals can be estimated 
with some confidence. This behaviour holds for all fields with ng < -3/2 and we refer to these as "coherent" fields. However, 
even for these fields the stress statistics become strongly dependant on the damping scale for k > kcoh- Since the damping scale 
is in general significantly smaller than the scales relevant for CMB analysis we can safely assume that the infra-red fields are 
"coherent" across the entire range of k. Only if small scales must be taken into account must more care be taken. In comparison 
to the case with rig - 0, for this field the cross-correlation between the scalar components obeys the same behaviour as the 
auto-correlations . 

To briefly summarise, then, we can consider the stresses of a white noise, decoherent field (and, by extension, fields with 
riB > -3/2) to be white noise on scales k < kc/lQQ. The stresses of a coherent field with ng - -5/2 (and, by extension, fields 
with hb < -3/2), become a pure power-law for k < kc, a far weaker condition that is easily satisfied for much of the history of 
the universe. 

As we demonstrate in the next section, only the largest scales are required to predict the CMB signals. On such scales the 
analytic solutions are significantly simplified. We expect ultra-violet fields to possess a well-defined Taylor series, with the 
zeroth-order term producing the white-noise regime. In contrast, we expect the infra-red fields to admit Laurent series with 
the leading-order contribution oc A:^"*+^^. In Appendix B] we present the large-scale expansions for all integer and half-integer 
riB e [-5/2, 3]. It is noteworthy that, for spectral indices ns > -3/2, the correlations to zeroth order are related to one another 
by 

Prr(k) 2Prr(k) ^Prrik) 5' 



as predicted for "causal" defect mechanisms in 116611 . (Note that the usages of the word "causal" differ between this paper and 
that work; here a "causal" field is one with hb >2 which can be produced by physical processes in a matter or radiation universe 
and retains the analyticity of the power spectrum.) 

In Figure |4] we present Prj and !P„ for hb e [-5/2,3]; the other auto-correlations follow the general behaviour of TrT- 
Consider Prr first, plotted in the left panel. It is immediately noticeable that fields with ng > -3/2 are decoherent, but below a 
"coherence scale" of A;coh ~ kc/lOO tend towards white noise. For riB < -3/2, kcoh - kc and both cases for which we possess 



' We employ the Numerical Recipes MISER routine. The complete set of solutions can be found in our notation in an online version of this paper found at the 



author's websitelhttp : //folk . uio . no/ibrownl 
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analytic solutions follow the expected scaling PAsik) ^ k^""^^ . For k > kcoh these fields are decoherent. At hb - -3/2 itself the 
stresses diverge logarithmically with k. 

The cross-correlation plotted in the right panel is of particular interest. In standard perturbation theory this cross-correlation 
is zero, and it is an interesting signal of a primordial magnetic field that it is in general non-vanishing. For ng > -3/2 the 
leading-order term is the gradient; the signal is thus negligible on the largest scales compared to the auto-correlations. Inter- 
estingly, however, the leading-order contribution to the cross-correlation is invariant with respect to hb and is, moreover, an 
anti-correlation. This feature is clear in the figure; while the scale at which the signal changes from an anti-correlation to a 
correlation varies slowly with ng the behaviour on large scales remains 

Prs = -^ABTTklk. (25) 

This behaviour is broken at n^ = -1 where the cross-correlation becomes positive-definite and on large scales behaves as 

Prs = ABTlklk. (26) 

The auto-correlations change their large-scale behaviour at n^ - -3/2, swapping from white-noise at n^ > -3/2 to a power law 
at hb < -3/2, and diverge logarithmically at hb - -3/2. In contrast, the cross-correlation is a power-law at both hb > -3/2 
and hb < -3/2, while it is white-noise at ng = -3/2. Furthermore, for ng < -3/2, the cross-correlation rapidly becomes of 
equivalent magnitude to the auto-correlations. If ng ~ -5/2, the existence of a non-vanishing scalar/traceless scalar correlation 
on the CMB is then an inevitable and characteristic consequence of a magnetic field. 
We now turn to two toy models with the power spectra 

rB(k) = Aek^ exp(-^k^)H( 1 - k) (27) 

with ^ > 20, corresponding to a "causal" magnetic field with a realistic damping tail - this can be compared for example with 
the power spectrum presented in lfT2ll - and 

Pg(k) = -Ab tan"' [^k^] k-^I^H{\ - k) (28) 

with ^ > 10^", which approximates a spectrum with hb - -5/2 for much of the range of k but is controlled on infra-red scales 
k < 10 '^. The damped causal spectrum will qualitatively resemble that for a magnetic field produced by physical processes 
in the post-inflationary universe - whether produced by induced currents, at a phase transition or through some other process. 
While the IR-controlled spectrum could be motivated by a magnetogenesis model producing fields on an extremely small scale 
before or during an inflationary epoch, more concretely this spectrum models a realisation of an inflationary field generated on 
a grid, with the infra-red damping reproducing the grid cut-off. 

Let us first consider the damped causal field, the stresses of which we expect to by similar to those of a power-law with 
riB = 2. Since kc remains a hard cut-off but the field is exponentially damped on a larger scale, we should see the gross features 
of the standard causal field move to smaller wavenumbers. We can solve equations (22p3 i analytically and the general forms 



are presented in the appendices and plotted in the left panel of Figure [5] The Taylor series around k = Q gives us the behaviour 
on large scales: 

!P,, = nABkl [To - ^k) , ^s = nABkl [-^k] , 
Tss = nABkl [lTo-4^k), Vy, = nABkl ( i| T^o - g^^) , . (29) 

Vrr^nABkl{^^^To-^k) 

Here erf (jc) is the error function and 

1 15 V27r erf( '^)^"'e^^ - 60^^^^ _ 39^9/2 _ 64^11/2 
^°-64 ^^VVi ■ (30) 



The stress spectra for a damped causal field retain the "causal" relationships given in equation (24i. Across the full range of 
k, the damped causal case qualitatively closely resembles the causal case; however, features seen in the causal case such as 
the flattening of the gradients iXk ■x. k^. are suppressed and there is a distinct loss of power on large scales. The behaviour 
of the scalar cross-correlation is particularly interesting. On extremely large scales the cross-correlation is an anti-correlation 
growing linearly with k, as for the standard causal case. This linear term is, however, strongly suppressed with respect to the 
quadratic term for k > 10"'^, above which the quadratic order dominates. Across almost the entire range oik, then, the cross- 
correlation is positive and growing quadratically. The level of suppression is dependant on ^ and in principle the behaviour of 
the cross-correlation on large scales would allow us to clearly distinguish between the causal and damped causal cases. 
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The IR-controlled field is conceptually simpler than the damped causal field, since the modifications are on very large scales 
rather than on scales comparable to k^. We would expect the stress power spectra for such a field to closely resemble the standard 
riB - -5/2 case closely above a certain wavenumber k^, but to decay on larger scales. Results from numerical integration are 
presented along with those for the standard case in the right panel of Figure l5l this spectrum is unfortunately not amenable 
to analytical integration. As expected, the small-scale results are almost exactly equivalent, although there is a consistent 
suppression of power for k < kc- On scales larger than a transition scale ks, however, the spectra flatten significantly compared 
to the undamped case. The most dramatic deviation from the previous behaviour is again in the scalar cross-correlation, which 
decays rapidly (!P„ ~ k^) as ^ — > 0. In any magnetogenesis scenario that produces an infra-red cut-off, then, we can expect large 
changes in the behaviour of the scalar cross-correlation. The auto-correlations are strongly damped but remain non-negligible. 
For this particular spectrum, the transition wavenumber is kg ~ 10""^, but this will be strongly dependant on the form of the 
damping. 

Comparison of the right panel of FigurelSlwith the realisations plotted in the right panel of Figure [3] shows a good qualitative 
agreement in the behaviour of the stresses. The damping effect is significantly more dramatic in the case of the realisations since 
the largest mode we can produce is only of order a//dim ~ <x ■ 10"^ where a > 4 to avoid aliasing; if required this could be 



modelled taking a smaller ^ in equation (28 i or enforcing a hard cut-ofF aXk x 10" . 

We close this section by briefly discussing the advantages and disadvantages of the three approaches we have taken. We have 
principally relied on analytical solutions. These have the obvious advantages of being exact and, in particular, allowing us to 
recover the large-scale behaviour directly from a Laurent or Taylor series. In an ideal situation, one would always use analytic 
solutions and we do so whenever possible. However, analysis is extremely limited. In this and previous studies solutions have 
only been found for magnetic power law spectra with integer and half-integer exponent n^, which lead to integer powers of k 
on large scales; here we have added to this a spectrum composed of an exponential and a power law. For the most realistic case 
of a field with iig x -5/2 this limits us to the two cases ng - -2 and rig - -5/2, which is not satisfactorary. Furthermore, 
since equations (22p3i have been derived employing Wick's theorem, the underlying magnetic fields must be Gaussian, which 



in most magnetogenesis scenarios is not to be guaranteed. 

The other tool we have extensively employed is the numerical integration of equations ( [22p3| ). These are relatively rapid, 
and even more so since - as we demonstrate in the next section - for purposes of CMB analysis we generally need only recover 
the region where k <^ kc- The approach that we have taken to the numerical integration helps control the various poles that can 
otherwise be problematic, and the agreement with the analytics is so good as to be indistinguishable. Furthermore, in a numerical 
integration we are not limited by the shape of the power spectrum, although for steeply-tilted spectra (see for example the case 
with rtg - -2.9 in the next section) we require an increasingly large number of samples in the integration to retain convergence. 
Numerical integration is therefore currently the most realistically flexible approach. Unfortunately, since it is based on equations 
(22p3 I we are still limited to considering Gaussian magnetic fields. 



The most general approach we have taken is to generate realisations of magnetic fields with a particular power spectrum and 
to construct and analyse the corresponding stress tensor The advantages are clear: this approach is entirely general. Neither the 
form of the power spectrum nor (in principle, although we have focussed on Gaussian fields) the statistics are at all constrained. 
However, the disadvantages are also clear. This approach is slow, involving numerous operations spanning the entire grid and 
frequent Fourier transforms. More seriously, results can depend strongly on the grid resolution. For ultra-violet power-law 
fields we can acheive reasonable accuracy, although the minimum wavenumber we can consider is limited by the grid resolution. 
However, for a power spectrum with a red tilt - with respect to the white noise case Vsik) - const. - there is a strong infra-red 
damping for all k < kc- This damping is lessened if we increase the resolution of the grid but this correspondingly slows the 
calculation further. It is possible to improve the results on large scales by employing a damping scale at k^ > /dim; due to the 
nature of the algorithms involved, power on smaller scales is strongly contaminated by aliasing from large scale modes, but the 
large modes themselves are less affected. Since it is precisely these scales that are most relevant for CMB studies, in cases where 
analytic or numerical solutions are lacking - when the fields are non-Gaussian, for example - this would be a crude method of 
finding the relevant behaviour of the spectra without running realisations on grids of size l^^ > 1 , 024, but it should be only be 
employed with caution. 

These issues only worsen when one considers instead the bispectra of the stress tensor Since the wavenumbers must form 
a closed triangle, on large scales we have only a small set of modes on the grid contributing to each configuration. While for 
power spectra we can recover reasonable results with (9(10) realisations, for the bispectra we are required to consider OilQr') 
realisations, as in B06. Worse, as a bispectrum is effectively a two-dimensional structure, the computation is significantly 
slower than for the effectively one-dimensional power spectrum. The need to employ such a large number of realisations in an 
unavoidably lengthy algorithm makes generating bispectra from realisations an intensive task. See B06 for further details, where 
we considered the one-dimensional colinear case; even for this simple case we required ~ 1, 500 realisations. 

Ideally, we would always employ analytical solutions. However, where this is not possible the most realistic approach remains 
numerical integration of equations (22p3 i. Although this currently restricts us to Gaussian magnetic fields this is not necessarily 



problematic where a concrete magnetogenesis mechanism predicting a specific statistical nature is lacking. In the remainder of 
this paper we will neglect statistical realisations and focus where possible on analytical solutions, although we employ numerical 
solutions in cases where these are lacking. 



10 

V. CMB ANGULAR POWER SPECTRA 

In this section we employ an approximation to the magnetised tensor transfer function to recover the CMB angular power 
spectrum induced by magnetised tensor perturbations in a general manner which will be directly applicable to the CMB angular 
bispectrum. This includes evaluations of the signals expected from the damped causal and the IR-controUed fields introduced 
in the previous section, but our focus is on the techniques and preperation for bispectrum evaluations rather than necessarily on 
the CMB angular power spectrum itself However, the same approach could be employed using numerical transfer functions (as 
in 147. .52.1 ) recovered from a magnetised Boltzmann code. While we consider only the temperature angular power spectrum 
arising from tensor perturbations, similar arguments apply to the polarisation and {TE) angular power spectra, and to scalar and 
vector modes. 

A. General Considerations 

The CMB temperature angular power spectrum can be found by integrating across the transfer functions, 

C, = ^ \nk) \^Ti{k, mf k^dk, (31) 

where in contrast to the normal definition we have taken the primordial power spectrum to be 

<^(k)r(p)) = nkxinfdik - p). (32) 

Here f (k) is a random variable of unit variance that characterises the primordial statistics of the perturbations. The transfer 
function for tensor modes is given fi3] l67l by 



where 



A^./(^,'7o)= J^^3^]^^^5< '(fc,,)^5^^-^d, (33) 



5^^^ = -2M^'+^<D(^' (34) 



is a source term, g is the visibility function and <1)^^' depends on the lower moments of the tensor Boltzmann hierarchy - see for 
example B06 for further details. 
Neglecting the photon and neutrino anisotropic stresses, the tensor perturbation evolves as 

a ■ T 167rG /^-^ 

h + 2-h + k^h^ — — T<^' (35) 

a a^ 

where an overdot denotes a derivative with respect to conformal time. After matter/radiation equality the source then rapidly 
becomes negligible. It can then be shown Il40ll5l1 that a good approximation for h in matter domination is 

h{k, ?7) * 327TGkT]lz,^ In ( — I ^^ (36) 

\Zeq/ kj] 

and then that, assuming instantaneous recombination, the tensor transfer function at the present day is approximately 



(^-2)! k^T^l 



^'M^o^-^ATTT^^^b^ (37) 



with Ji(x) a Bessel function of the first kind and D a constant. Employing this expression assumes that we are working on large 
scales and it should therefore only be trusted for relatively low multipoles. 

At the present epoch the transfer function is extremely tilted to low k and we are restricted tok < 10"^- 10"'*Mpc"' . In contrast, 
the damping scale relevant for tensor modes is kdrj^q) which is of the order of a few to a few hundred inverse megaparsecs. The 
transfer function then restricts us to extremely large scales: for inflationary fields we are well into the power-law regime, while 
for ultra-violet fields the sources are close to white-noise. The damped causal field also produces a white noise source on such 
scales. In general, we expect the majority of field configurations to produce source spectra that on large scales act as 

'PTT(k)^pA^\ (38) 
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for some amplitude f * and pivot scale k-i,. For ultra-violet fields (including the damped causal case) we have a - 0, while for 
infra-red fields a = 2nB + 3. In the case of fields for which we possess an analytic solution we can further identify P-t/k^l from 
the Laurent series (see Appendix [B|. In a more general case, given sample points from a numerical integration or the results 
of an ensemble of realisations, we can recover the scaling a and amplitude P^, around a pivot wavenumber A:*, which we can 
employ in the CMB integration. The structure of the power spectrum on small scales is practically irrelevant so long as equation 
( [38] l holds while the transfer function is non-negligible, which we expect to be true for all magnetic power spectra that do not 
possess features in the deep infra-red. 



In principle the integral in equation ( 3 1 1 runs only to ^/fef (?7eq) - 1 since above this scale the stress spectra become decoherent. 



regardless of hb- However, since the transfer function is so steeply-tilted, the error introduced by integrating to infinity is 
negligible for any realistic 'PrTik). Employing the approximation (38 1 then gives 

m nD^ P^l(l + 2)\ f" „ . , 
This is a standard integral and the solution is 

rn'D'p. i r(^)z(/ + 2)!r(^^) 

This is true for all systems in which the approximate form for Aj-; holds, and for which 'Prrik) obeys a simple power-law while 
AL is non-negligible. Since in evaluating the numerical results in section IV we set Ag = I and worked with a wavenumber in 



units of kc, then if we recover a value P* at a pivot scale A;* = k^kc we have that 



P^Agkl 'P„={ABkl-''P^)i:^] . (41) 



We should emphasise that equation (40 1 has been derived based on the form oi An, which equally applies when considering 
the magnetic bispectra ll68l . The magnetic bispectrum is a 3D quantity dependent on wavenumbers {k, p, q] which form a closed 
triangle. We can parameterise this with the variables {k, r, (p) where r - p/k and (p is the angle between k and p. It has been 
shown that with r - 1, for = (the so-called "colinear" case P3]|59ll6T| ), and for both (p = 27r/3 (the equilateral case) and 
(p ^ JT (the degenerate, or squeezed case) ll60ll6Ti the bispectra for «g < -1 follow scaling laws in kfoik< kcoh- In forthcoming 
work we further demonstrate that the coherence scale for = is A:coh ~ ^c/2; for greater values of <p we expect A:coh ^ kc/2. The 
nature of Ati then implies that we can employ the same reasoning as we have for the power spectra and focus only on extremely 
large scales, sampling only as many points as necessary to reconstruct the behaviour of the intrinsic bispectrum SxTTik)- While 
the benefits of doing so for the angular power spectra are arguable, the saving in computational time and complexity for the 
bispectra will be extremely significant. We explore these issues thoroughly in a follow-up paper 

B. Results 

Consider first ultra-violet fields, with riB > -3/2. In this case the power spectra are white noise on the largest scales and 

a = 0, P^^ TxABklfiriB) (42) 



where /(ng) is the scaled amplitude of the zeroth-order term in the Taylor series. The CMB angular power spectrum (40i is 
therefore 

^ _ 16;r^ ^2„2 lkc^\' X^'firiB) Kl + 2)! {21 + 1)!! 

^'- 3 ^""''Uoj r(^) (/-2)!(2/ + 9)!!- ^^^^ 

Here a\\ - a{a - 2){a - 4) ... is the double factoral function. In the limit of large / the angular power spectrum tends towards 

/(Z-Hl)C, oc/3 (44) 

in agreement with earlier estimates BOl . In particular, consider the causal case with ub - 2 and the flat case with ub - 0. For 
hb - 2, /(ng) = 16/15 and the CMB angular power spectrum is 

l,0247r5/2 ^lkcX\' .l{l + 2)\{2l+l)\\ 

Ci^- D^B\\^\ A^- —. (45) 

135 Wo/ (/-2)! (2Z + 9)!! 



For hb = 0, fins) - 1 12/45 and so 
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135 \i]qI (/-2)! (2Z + 9)!! 2A^ \kc(ri^^,nB ^2)j 

Since the damping scale for ng = is significantly greater than that for ng - 2 BOllSTI . all other parameters being equal the 
signal from a flat field will be larger than that for the causal field. 

The situation is different for the infra-red fields as a becomes ng-dependent. For hb < -3/2, 



3,58471^/2 2 /MV /(' + 2)! (2/+ 1)!! _ 7 /^c(?7eq,nB = 0)\^ 



a^2nB + 3, P*^„« = nABk^,-"finB) (47) 



and ultimately 



c = —d'b' f^f -^"^/("ii) r(^-"^) ^a + 2)! r(/ + 3 + «g) 
' 4 ''Uo/ (;fc,,,o)2«''+3r(^)r(i-«B) (/-2)!r(z + 4-«B)' 

This holds for all power-law spectra with hb < -3/2 but since ng is not necessarily an integer or half-integer the dependence on 
/ is not necessarily straightforward. 

There are two fields in this regime for which we have exact solutions. If n^ - -2 then /(n^) = 9jP' lA and the CMB angular 
power spectrum takes the simple form 

Ci = D^B]—^ ^ (49) 

128 -* T]l (l + 5)(l + 4)(l + 3) 



while for hb - -5/2, we have /(ng) = 3, 008/75 and 

770, 048 
l,125r(l/4) '' 770 (^-2)! (2Z-H11)!! 



C, = ^^^^^^.^'Bl'^^Sil^^^. (50) 



In the limit of large-/ these become 



'^'^i>^"^^'/;"'=:5/2' (51) 



again in agreement with previous results BOll . Note that these results have been found with relative ease from a flexible approach. 
We also have an exact solution for the damped causal field. The amplitude of the spectrum is given by 

'" - 8.2(/;,2p/2 (3 V^ erf {kA^ . A^y^) - (6 . Ak^ (^ . A^)) U^ . ^^y /2^-,^(,..^)) (5^) 



which for reasonable parameters (kcA > 5) is rapidly dominated by the error function giving 

\5/2 



87T^''-{^ + A^f 



Ab * '— —Bi (53) 



From the stress spectrum in equation ( 29 1, on large scales we have a - and 



^3 7 15V2^^Vf erf (V2I)- 4(15^^/2 + 20^''/2 + 16^1 1/2) 



^* - -^'^''WO .13/2,2, ■ (54) 



The CMB angular power spectrum is then 



r 28.rV2 ^,,5/2 /^, -|^ /(/ + 2)! (2/ + 1)!! f 15 V2^^V^ erf ( V2|) - 4(l5^^/2 + 20^^/^ ^ 16^"/2) -| 
^'"^35-^^-^!^^'^'' y (/-2)!(2/ + 9)!![ e^JV^ J- (^^^ 

In the limit ^ — > this reduces to the angular power spectrum for an undamped causal field ( [45] l, while for large ^ it becomes 

C. . ^ ^D^B] (kPll^ (^[±i)!l . 1 ( ^(^i^^^^^P^)^ CKundamped). (56) 

9 ^ ''[rjoj (1-2)1(21 + 9)11 f^^ \^/t,. (undamped)/ '^ t' ^ ' ^ 
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As should be expected since, on the scales where A^/ is non-negligible the stress spectrum for the damped case is similar to that 
for the undamped case, the CMB signals are similar between the two cases. 

Consider the "transitionary" case with tig - -3/2. For this field, the stress spectrum on large scales diverges logarithmically 
and we cannot necessarily employ a simple scaling approximation for Vrr- The amplitude of the power spectrum is given by 
equation (J7]i. Employing the approximation in equation (B12i and the general expression in equation (31 1 gives, after some 
manipulation. 



C, 



135r(3/4) 



©^B^ 



'70 



ln(fe,77o)-*P U+x +ln(2)- 



1 /(O {21+ 1)!!U(/ + 2)!(2Z+ 1)!! 



6 2Z+ 1 (2/ + 9)!!/ (/-2)! (2/ + 9) 



(57) 



where /(/) = 1607^ + 2,576/"* + 15, 152Z^ + 48,618/ + 20,529. Of the terms summed in the brackets, the most significant 
contribution come from the first two; for sensible parameters we always expect have kcrjQ » 1, while the term dependent on the 
polynomial /(/) decreases monotonically and is only marginally significant even at / x; 2. We can therefore approximate 



Ci 



\J91t? 
135r(3/4) 



D^BJlA 



..^)-(,„(,,,„)_^(,4 



/(/ + 2)!(2/+l)!! 
(1-2)] (2/ + 9)!! 



(58) 



although in Figure l6] we retain the full generality. In these expressions ^'(x) = d(ln(r(x)))/dx is the digamma function, which 
is potentially significant for larger / > 1000 but negligible for larger angular scales. Note also that on the scales on which *i'(x) 
might dominate the approximations for the transfer function will no longer be valid. The impact from this field then closely 
resembles that for ultra-violet fields with hb > -3/2, although the similarity grows less for increasing /. This is to be expected 
since the divergence at large scales in Vjjik) is only logarithmic and to some extent could be approximated by white noise. 
The modifications to the CMB angular power spectrum at large scales tend to decrease the /-dependance of the signal; this is 
consistent with the general flattening of /(/ -i- 1)C/ as ub reduces and grows more negative. 



The magnetic spectrum for an "inflationary" field controlled on large scales, equation ( 28 1, was chosen to act as an ultra-violet 
field on large scales but to mimic a field with hb - -5/2 on smaller scales. This induces a stress spectrum split into two regimes, 
characterised by a change in the gradient on large scales. If the gradient is a\ on scales k < ks and ai on scales k > kg, then 
assuming ^* <ks. 



^ttW- ,a, i k"'-"'-k"^- ,ke(ks,k,] ■ 



(59) 



From the numerical results shown in Figure [5] for the IR-controlled field we can take a\ - Q and q;2 - -2, the deep infra-red 
acting as white noise while smaller scales obey the same law a2 = 2«b -i- 3 as the infra-red fields. Substituting this into the 



general expression for the magnetised CMB in equation ( 3 1 1 and solving the integral yields 



n .P^ .1(1 + 2)1 



4 tjI '(1-2) 



512 (21- 1)! 



1 



1571 (2/ +11)!! 32[(/ + 3)!]2 



(ff" 



1 



(21 + 3) 



r. 



2/+7,/+4,/+| 



(-<) 



1 



T. 



1+1,1+1: 



(2/ +!)(/- 1) 21+7,1+4,1+ ^_ 



A-x]) 



(60) 



where 'F^^f^''(x) is a hypergeometric function and x, - ksTjo. Corrections from the hypergeometric functions are restricted to 
very low /, and unless ksTjo > 1 are negligible. Since we expect ksfjo <*: 1 by construction we have 



15 tjI (l-2)\ (2/+ 11)!! 



(61) 



Even though the stress spectrum decays on large scales, where the transfer function is most significant, this is indistinguishable 
in form from the imprint of the straight power law with n^ = -5/2. Writing the magnetic power spectrum as 



'PB(k) = ABk-^l^ tan"' [^k^)H(kc - k) 



(62) 



then the normalisation (J5]l has an unwieldy exact solution dependant on a number of hypergeometric functions. However, we are 
assuming that ^ is very large implying that the hypergeometric terms are negligible and the field can be normalised to 
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P, 
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^expected 


-2.9 


9.233 X 10" 


-2.80 ±0.04 


-2.80 


-2.75 


4.552 X 10'" 


-2.50 ±0.02 


-2.5 


-2.5 


5.024 X 10** 


-2.00 ± 0.02 


_2 





7.812 


(-2.47 ± 7.40) X 10-"* 





2 


3.345 


(-7.76 ± 8.38) X 10-* 





2 (UV controlled) 


9.627 X 10-5 


(-0.49 ± 1.36) X 10-3 


~0 



TABLE I: Stress spectra recovered for sample power-law magnetic spectra. 



From the numerical integration, the amplitude P^, evaluated at a pivot of A:* 
One can then find that the CMB angular power spectrum is 



k^/k, ^ 7x10-*^ is P* = P^Aekl ^ (4.5xlO^)AB/t^. 



C, 



: 46P^klD^BJ 



A^'^kl 1(1 + 2)1 (21- 

m (i - 2)! (27T 



1)!! //t,.(IR controlled) \^ 



where in the second step we have used that ks ~ fe^/lOOO. Since the damping scales are approximately of the same order 
of magnitude, the signal arising from such an infra-red controlled magnetic field is smaller than, but otherwise practically 
indistinguishable from, the undamped case. 

The methods we have employed for the IR-controlled field demonstrate a general approach to the CMB angular power spectra 
applicable for any magnetic field for which we can numerically generate - or produce with statistical realisations - a stress 
spectrum. In the simplest cases we need only recover P* at a pivot fc*, find a, and then employ these in equation (WOb. TableH^ 
shows the values of P* and a and the l-cr errors on a, for a range of spectra, with ki, - (5 x 10-'*)A;c(77eq),'* We recover a from 
averaging the gradient of the numerical spectra across much of the coherent range - typically from ~ 20 samples. Figurel6]shows 
the resulting approximations. Note that before recovering P* and a we have smoothed the numerical integration for «b = -2.9 
to remove some residual noise arising from the limited sample volume. Since the scaling regime is clear, and only the accuracy 
with which a is recovered depends on the smoothing, there is no error introduced in doing so. The saving in integration time 
when we need only consider a limited range in k is significant even at the two-point level; when we consider the bispectrum in 
future study this will become even more important. 

By extension, we can model more complex spectra, as with the IR-controlled field, with an arbitrary number of distinct 
regimes, each characterised by a particular value for a. However, in most cases we would expect to require only one or two 
regimes. (An alternative approach employing general fits for ranges of ng was recently taken in |27 1.) Similar arguments will 
equally apply to the vector and scalar perturbations and to the polarisation power spectra, although the form of the transfer 
functions will be different. More directly, similar arguments apply for the tensor CMB angular bispectrum, which can be 
recovered employing the same transfer function. 

To close this brief consideration of the magnetised CMB we consider two power-law fields, with ng - -2.9 and ng = -2.75. 
For these fields we can employ equation (40 1 with the values of P* and a from Table Q Employing the central value of a and 



neglecting the error we can wrap the infra-red power spectra onto the microwave background with equation (40 1 to find 
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Analytic 



hb = -2.5, 
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(66) 



** Note that the physical scale at which we are sampling is not constant; we are choosing a pivot point in units of kc = kcirj^i 
therefore employing a pivot on much larger scales than we are for lower ng. 



, «s). For large values of «b, we are 
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_ n^l^ YQAm , p.B^^2„2M'^^^f f ° r(/+ 1/10) 

' 4r(i/20)r(39/io)^ * * '' ^\ m I r(z + 69/io) 

, ,//i'/2;tf\'^"' r(/+i/io) 

^ 90.9AD'b\\ ^ ; ' ' ub^ -2.9. (67) 

'*\ 770 / r(/ + 69/10) 

The agreement between the numerical and analytic estimates for ub - -2.5 is extremely good and confirms the approach, while 
the /-dependence for the other two cases is in agreement with previous estimates within their range of validity (for example 
||40l|55]|57l). 

In Figure l6] we plot the large-scale CMB spectra for the magnetic fields we have considered in this paper normalised to 
1(1+ 1)C/|;=2 = 1 to emphasise the /-dependences. In this plot we employ the full result for hb = -3/2 field rather than 
neglecting the *!'(/ -1- 1/2) and /(/) terms. We would emphasise that it is not our intention to derive exact CMB predictions - 
that has been done in, for example, 12611271 52, 55, 57 1 employing the results from magnetised Boltzmann codes. It has rather 
been our intention to demonstrate a rapid approach to the large-scale magnetised CMB that will prove extremely useful when 
considering the magnetised bispectra, as it allows us to focus purely on the k <K kc region. However, we have here presented the 
impact on the microwave background of two interesting fields, one a causally-generated field damped in the ultra-violet across 
a wide range of k, and one an acausal field damped in the deep infra-red. In both cases, for the reasonable parameters we have 
chosen, the impacts are indistinguishable from a standard causal field and a standard inflationary field respectively. 

VI. CONCLUSIONS 

In this paper we have considered in detail the statistics of the stress tensor of a tangled primordial magnetic field at the one- 
and two-point levels, employing a mixture of analytical techniques, numerical integration, and statistical realisations. Assuming 
Gaussian statistics and a power-law spectrum for the magnetic field, we are restricted to the region n^ > -3; a field with n^ - -3 
would be exactly scale-invariant. Fields with ub > 2 can be produced by causal processes. We have examined the one-point 
moments, extending the results of BC05 and B06 to a wide range of spectral indices and a much higher grid-resolution. We 
have verified that the isotropic pressure of the magnetic field is a x^ field to a very high degree, as has been argued before. We 
also considered the anisotropic pressure, which runs with spectral index. In particular, for indices n^ ^ -3/2 the statistics of 
the anisotropic pressure are constant, with a skewness y\ x -0.24 and a kurtosis 72 ~ l-l- However, in the region riB < -3/2 
both the skewness and the kurtosis increase as hb — > -3. There is a strong grid dependence for such steeply-tilted spectra which 
requires us to operate on grids of side-length /dim = 512, which arises from the power being piled onto a small number of modes 
on the finite grid. This infra-red divergence is inevitable - and re-emerges in the study of the 2-point moments - but we can 
nevertheless predict that yi < 1 and 72 ^ 2 at n^ = -3. Perhaps more importantly, we have demonstrated that at the one-point 
level the field with hb - -3/2 separates the parameter space into two regions in which the anisotropic stress behaves quite 
differently. 

Turning to the two-point moments, our chief focus has been to construct a formalism which can be equally applied to both 
analytic and to numeric evaluation of the stress spectra - and which, more importantly, will be readily adaptable to the calculation 
of magnetic bispectra. We then applied this to the two-point moments of a power-law magnetic field, finding analytic solutions 
existing for all integer and half-integer hb e [-5/2, 3]. While some of these solutions were previously found in |57 1, others were 
not. As has been known for some time, the parameter space is split in two, again at hb = -3/2. For hb > -3/2 the integrals are 
readily solved, and across the range of k, the stress spectra are technically "decoherent", in that the amplitude and the gradient 
at each k depend explicitly on the damping scale, which is time-dependent. On the very largest scales, for k < kcoh ~ ^c/100 the 
cross-correlation of the isotropic and anisotropic pressures vanishes (and is an anti-correlation on intermediate scales, except at 
hb - -1), while the other correlations tend towards white noise. The ratio between the various correlations is constant across hb, 
and is equivalent to the "causal" ratios found in |66|. Conversely, for hb < -3/2, there is a large region in which the statistics are 
"coherent", for wavenumbers up to A:coh ~ ^c- On these scales the power spectra scale as Pik) oc ^2«fl+3 ^^ jj^^ been appreciated 
for some time; it is worth noting that it is only for integer and half-integer hb that this results in integer powers of k. Interestingly, 
the cross-correlation between the scalar pressures is of the same order of magnitude as the auto-correlations. Above ^coh the 
statistics remain entirely decoherent. In particular, we presented a solution for the case with hb - -2. Since we cannot find a 
solution at hb - -3, there are only two solutions in this "coherent" or "infra-red" region available, hb = -5/2 and hb - -2. 
Since these values remain in the centre of the allowed constraints, it is important that we possess the solutions. Finally, the field 
with hb - -3/2 behaves quite differently. In this case the cross-correlation of the scalar pressures is a positive correlation, and 
white noise on large scales, while the auto-correlations all diverge logarithmically with k. Numerical integrations employing 
Monte-Carlo techniques confirmed the analytic solutions to an extremely good accuracy. 

From our analysis, we then not only have an approach to analytical and numerical integrations that can be extended with ease 
to bispectra, but we have also clearly identified the important regions of the spectra solutions, broadly governed by a coherence 
scale ^coh- 
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We then considered two non-power law magnetic power spectra - the first, a "damped causal" field has a spectral index 
«B = 2 on large scales but an extended damping tail on smaller scales, and resembles a field produced by a causal mechanism 
such as the nonlinear vorticity of the electron/proton plasmas; while the second, an "IR-controlled" field, has a spectral index 
hb - -5/2 on relatively small scales but decays on large scales, and while it could in principle be produced by a carefully-chosen 
inflationary scenario is more concretely employed to model the infra-red cut-off of a sample realisation. Analytic solutions exist 
for the damped causal field we chose. The damped causal stress spectra resemble the standard causal case with the gross features 
moved to a smaller wavenumber corresponding to the peak of the magnetic power spectrum. The IR controlled spectra, in 
contrast, tend towards white-noise for k < ks where ks is related to the scale at which large-scale magnetic power begins to 
decay, but closely resemble the standard infra-red case on larger scales, albeit with a suppression of power for all k < kcoh- For 
both of these scenarios, though, the cross-correlation significantly differs from the undamped cases, which might in principle 
allow us to distinguish between them. 

We also considered the power-law spectra employing statistical realisations, running on grids of side-length l^i^ = 256. While 
the results in the ultra-violet regime with hb > -3/2 are in extremely good agreement with the analytical solutions, the limitations 
of the approach are quickly apparent: we are limited to a minimum wavenumber k/kc ~ 4//dim ~ 10^^-10"^. While this is close 
to the white-noise regime, one would not necessarily wish to naively rely on the output at such a large scale, particularly as the 
error at the minimum wavenumber is naturally rather large. When considering hb < -3/2 the limitations of the realisations are 
even clearer; there is an extremely strong impact from the infra-red damping. While working on increasingly large grids eases 
the problem this rapidly becomes unwieldy and unrealistic. While realisations remain in principle the most flexible approach to 
the statistics of early universe magnetic fields, not least as they can easily take into account magnetic fields with non-Gaussian 
statistics, they should be treated with caution. 

Finally, we outlined an approach to the CMB which will be of use in the study of the CMB angular bispectrum arising from 
magnetic perturbations. Working with the tensor perturbations sourced by a magnetic field we considered the temperature auto- 
correlation. Using that the transfer function is steeply tilted to large scales, we argued that generally only the regions k < kcoh 
contribute to the CMB integrations. The stress spectrum can then with fair generality be written as a simple power law, where 
the amplitude P* at a convenient pivot scale A;* can be recovered from an analytic solution, from a numerical integration, or from 
a statistical realisation. Assuming that the approximate form for the tensor transfer function is accurate, exact general solutions 
for the CMB signals can then be found, and specialised to particular power spectra with a suitable choice of P* and a. Where 
this is not true, an approach similar to that which we took for the IR-controlled field can be applied: the stress spectrum can 
be split into two or more regions, each characterised by a different a,. The solutions that result are unwieldy combinations of 
hypergeometric functions, but in most cases are likely to reduce to simpler forms for realistic parameters. 

In addition to reproducing the previously-known forms for the CMB angular power spectrum from magnetised tensor modes 
and power-law magnetic fields, we considered the damped causal and the IR-controlled fields. For the parameters we chose 
in this paper, which seem reasonable, the signals turn out to be indistinguishable from those from the standard causal and 
inflationary cases respectively. However, we have presented the expressions for a wider variety of parameters. We have also 
evaluated the amplitudes of such fields when smoothed on a comoving scale k^. 

This approach to the CMB entirely separates the form of the stress spectrum from the transfer functions in a transparent 
manner. Furthermore, the arguments we presented do not rely on the specific analytic forms of the stress spectra; rather, they 
rely on the form of the magnetised transfer functions. Naturally, these do not differ between the angular power spectra and the 
angular bispectra. The bispectrum can be described as with the coordinate system {k, r, 0) and foliated into planes of {k, r); in 
each of the ID lines of constant-r through these planes it appears 1*451 l59]j6Tl that the bispectra on large scales behave as either 
white-noise in k or with a power law in k. Our approach to the CMB angular power spectra, therefore, can equally be applied to 
the magnetised CMB angular bispectra, which are otherwise extremely challenging to consider. This will be the focus of future 
study. 
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Appendix A: Analytic Power Spectra 

Analytic solutions exist at least for all integer and half-integer ng > -3. For integer ng they take the form of simple polyno- 
mials, while for half-integers they include logarithmic and inverse trigonometric terms. Here we present these for hb e [-5/2, 3] 
that have not been previously presented in [i57J . covering the magnetic fields range of most interest. In the following equations, 
we define 
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with equivalent expressions for the other correlations. In the following subsections we tabulate the non-vanishing J[„,{k), S(k) 
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2. A Damped Causal Field 



For the power spectrum defined in equation ( 27 1 analytic solutions can also be found. Here we define 
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^'i = C^i = i^(4(5-2A2(^))^2+4(io-7A2(^))^ + 31-10A2(^)), 

^o' - Cf,' = -^(1 - 12^)A2(^) - I V2;^.-^'^^/2£/(fe), (^^^^ 

y[f = Cf = -3|3(4(2+A2(fe))^ + 8 + 3A2(fe)), 

^f = Cf = -^,-f*^(^+ V2^,f(^^+4)/2£^(fe)), 

^f = <^f = 6|i(l-2A2(^)), 
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^- = C"; = -|i(l-A2«), 

^^3 = CI^s = ^ ((1 + 1 lA2(m + 2(1 + 2A2(fc))) , 
" A^i \ (A14) 

^w ^ qv ^ _^(2^ + 2+A2(fe)), 



^"5 = ^""5 = |^(l-A2(fe)), 

^"4 = ^"4 = |^^2(fe), 

^"3 = C':^3 = ^(2(l-4A2(fe))^ + 4-7A2W), 

^"2 = C':; = ^(7+4^)A2(fe), 



2 ~ ^^-2 ~ 2f* 

^"1 = C':^i = lli[A{\-AA2{k))e + ^{l-5A2{m-{\ + \0A2{k))), 

Jl^ = C- = ^(l+4^)-^e-f*^/2£^(fe), 

^r = ^7 = |T(2(l-2A2(fc))^ + 2+A2(fe)), 

y[^ = C^ = -^A2(fe), 

^" = C" = 3|i(l-4A2(fe)), 

3. Hb = 'ill 



(A15) 



yil'; 


- 


c- 


= 


^{^-^-"')^ 


mi^ 


- 


C5^ 


= 


1,226-1/2 
1,617" 


jiY 


= 


C7 


= 


4 /, 1,811 -l/2^ 

21 V ^ 308 " y 


jiY 


- 


q- 


= 


2,693-1/2 
4,312" 


:hy 


= 


CY 


= 


1 (-[ 3,449 -l/2^ 
9 l,^ 8,624" ; ' 


^7 


= 


CY 


= 


1,927 1/2 
44,352" 



^["3 


= 


C-3 


- 


i^A'-x-'")' 


-^-2 


- 


C-2 


- 


16 -1/2 
1,155" 


y["l 


= 


C", 


= 


244 /i 312 -l/2^ 
1,617 \^ 305" ) 


^0 


= 


C5^ 


- 


596 
8,085 ' 


y[7 


= 


CY 


= 


A/l _ 818 -l/2\ 
21 V 385" 7' 


^2' 


= 


C2 


= 


737-1/2 
980 " 


yl^^ 


= 


CY 


= 


1 /, 62,903 -l/2^ 

18 V ^ 10,780" ; 



(A16) 



aiTT _ ^TT _ 8,797 -1/2 airs _ ^ts _ 1,079 -1/2 

"^5 ~ '-'5 ~ 88,704" ' "^'4 - '-'4 ~ 110,880" 

^r = -5§6(^"'« + ^)' ^? - C-' = 25f50^-'/\ 

Q ~ ~2i>88'^°^ •^^)' 
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with 



Jl^l 


= 


c'l 


Jlil 


- 


c'_l 


J?[^f 


= 


c'l 


^-l 


= 


c- 


J?[^f 


= 


C'J, 


<ziSS 

■^0 


= 


'-'0 


J?[f 


— 


Cf 


J?[^^ 


— 


Cf 


^f 


- 


cl' 


^f 


— 


cl' 


y[f 


- 


cl^ 




^f = 



c. 



ss 



1,536 
80,465 

768 
80,465' 
128 
1,155 



[l-K-'^^)., 



-1/2 



(1-^--'^^)^ 



13,088 -2 
'241,395'' 



(>- 



223,942 -1/2 



J76 

147 V" 218,405'' 

2,003,102 -1/2 

1,689,765'*' 



Y 



21 V 
_ 275,949 

409,640 ' 
_ _1_ 

36 
_ 19,583 



2,053,913 -1/2' 
321,860 '^ 
-1/2 



:(' 



4,213,440' 



1,182,337 - 
204,820 '' 
-1/2 



104,837 -1/2 
8,426,880" 



155 l^T^-i 
2,688 



(cos 



(k) + n) 



59,136 



COS {k). 



)^ 



\I2\ 



rnw _ r>w _ 



JXry 



^7 
J{7 



2,048 
"241,395 



(l-^-'^'), 



^'-4 




^-A 






jr^ 


= 


C^^ 


- 


3 


-^-2 


= 


C'-i 


- 




jr, 


= 


c"; 


- 


296 
1,617 




= 


'-'0 


- 




y[7 


= 


CY 


- 


2 
21 



64 
3,465 



(> 



(1 
(' 



a 






- ^4 



^- = 



1,024 
'241,395 

- ?^jf 
209" 

6,304 

724,185' 

114,742 - 

115,995'' 

1,207,628 

1,689,765* 

866,661 - 

80,465 '' 

38,609 

55,860* 

249,103 

"1,228,920' 

403 

33,440' 

12,519 



l/2~\^ 
l/2\ 



1/2 



(A17) 



65 



29,568 ' 



a^'^ 



29,568'- 



468,160' 



.-1/2 






-2 









JCi 



CI 
C 



-3 



^rr ^ fjTT ^ 






1,024 
241,395 
512 ,,-1/2 



(1-^-'^'), 



241,395 
128 
3,465 



(1 



206^-1/2 
209 



13,568 -1/2 
724,185" 

104 /i _ 37,724 
1,617 V 40,755' 
2,161,256 -1/2 



.■'i 



1,689,765 



K 



_ /^7T _ 4 /1 1 

- '-'1 - 21 l^ 



,324,971 , 



rnTT /^TT _ 



c? 






857,909 
307,230* 

1^(1- 
49,501 



80,465 
1/2 



-l/2^ 

-l/2\ 



(A18) 






305 



C7 = - 



14,784 
305 



1,053,360' 
35,737 -1/2 
300,960" 
-1 
COS 



1,521,059 -l/2'> 
102,410 " / 

-1/2 



(cos {k) + 7n , 



14,784" 



'(k), 



(f>TT _ 503 _ 
■" ~ 4,851 



<gSS 



5,324,279 
20,277,180 



S"' ^ 



10,752^' 

155 
118,272"' 

3,564,031 



gv 



305_ 

10,138,590 29,568 



2,879 

48,510 5,376 
_ 495,794 
~ 5,069,295 

n. 



59,136' 



(A19) 



4. «B = l/2 



W\ = 



-^0 



rnTT _ _ 103 4-1/2 
■-^'2 ~ 360 " ' 



(1-.'/^), c- = 1(1 

^1/2 /-TT _ _277. 

" ' '-'o ~ 225" 

= -1(1- li-'^^)' cr = -1(1- 



104 

225 

123^^1/2 
225" ' 



+ K 



-1/2 



277^-1/2 
225" 



)^- 



289^-1/2 



C7 = 



1,800 



360 

1/2 



)' 



(A20) 



cr = (i + l-"^')' 



2,048 -1/2 
459 



K 



22 



with 



rnrs _ 
rnrs _ 



-2 



J?[" 



(JITS _ 

^7 = 



(i-/^-'/'), 



32 
195 

195 « 
148 /i _ 496 „-l 
225 l^ 481 



(1-if-"^^)- 



<2IT5 _ 



.932_ -1/2 
2,925" 



l,^ 2,925" }■ 



3,878 -1/2 

2,925 
7J47 1/2 
5,850 



K 



-Hi _ ^23_ -l/2^ 
2V 1,170" 7' 






^-2 



'40 ' 



C"l - 



q^ = 4(1 + 



^-1/2 



32 
195 

195 « 

148 

225 
J32_ -1/2 
2,925" 



(>- 



MZI;,-1/2A 

2,925" 7' 



(A21) 



.7447-1/2 
5,850" 



1,170' 



-1/2 



)' 



(55S-1(^) + 2i2^-l/2)^4^ c^. ^ |. (~-l(^^ ., 2^^-1/2) 



■ss 
•ss 
•ss 
ss 
ss 



aiSS 
^0 









SS 

ss 
ss 
ss 

4 



3,315(1 K jk 

-1/2 



256 

3,315^ 

128 

195 



(> 



« •'-), 



35^-1/2 
34 



J52^^-l/2 



1,105 
356 /'i 20,614 -1/2 



(> 

123 
'25 



19,669' 



107,123 -1/2 
49,725 " 



^ 79,560" 



43,057 -1/2 
397,800" 



- 1(1- M05 .- 

~ 4 \^ 7,956" 



1/2 



)' 

'^^), 

)> 



^(£5S-'(fe) + 5i«>?. 



160 



17,901 ' 



-1/2 



r^ss _ 
'-'-4 - 

•^0 - 






m{^^ ><-"')' 



256 ^-1/2 



3,315 
128 
"195 

352_^-l/2 



(i + i--'^^)^ 



1,105 
356 

225 



(1 



20614 -l/2\ 
19,669" ; 



Cf --l{l- 



107,123 -1/2 

49,725 " 

72,071 
79,560* 



(A22) 



-1/2 



43,057 -1/2 



)■ 



397,800 



K 



4V ^ 7,956" 7' 



), cf = -4(£5s-w + ?i^. 



17,901 ' 



l/2\ 



<7JVY _ 

<7IVT/ — 

rztVV _ 

A) - 

<7IVV _ 






2,048 
■ 29,835 
1,024 -1/2 
29,835" 



(1-^-'^'), 



64 
585 



(>- 



51' 



1,504 1/2 
29,835" 



(>- 



104 

225 

174,446 

149,175" 

32,641 - 



8,414 
8,619 
1/2 



'^^), 
^-'^'), 



/c 



198,900 
1.067 ,,-1/2 
26,520" 



41,737 -1/2' 

19,890 

1/2 



~ 240 V 



COS (k) ■ 



2,048 ,,-12' 



663 



K 



C!'5 = 



C!: 



2,048 
'29,835 



1,024 -1/2 
29,835" 

585 (1 + 51'^ 



(1+^-'/'), 



1,504 -2 



C7 = - 

cr = - 



29,835 
104 
225 

174,446 . -1/2 



(> 



^ 8,619" 



l/2~\^ 
l/2\ 



(A23) 



149,175 



K 



(1 



41,737 l/2'l 
19,890" ) 



32,641 1/2 
198,900" 



2\ /^w _ 

j' <-4 - 



1.067 ,,-1/2 
26,520" 



^"5 = 

<ZI7T _ 



-3 



^1 



^0 - 
A - 

O?^ — 

^7 = 



1,024 
29,835 
512 ,,-1/2 



(1-^-'/'), 



29,835 

128 

585 



(l-i--'^^)> 



J56. -2 
2,295" 



251,468 ,-1/2 
149,175" 

4/, 34:031 -l/2\ 
5 l>^ 6,630 " ) ' 
iMZ3.-l/2 
7,650 " 
1 /-i _ 3,977-1/2 



6,630"^ 
'-1/ 



)' 



C!_'5 = 



c 



-2 



- 88 /i 7,828 -l/2'\ ^rr _ 

~ 225 V 7,293" 7' '-'-1 ~ 



c? 



(1 



UHl/i -l/2\ 

29,835 \^ ^" ; 
512 -1/2 
29,835" 
128/1 , 50 



585 



(> 



+ f^K 



-1/2 



_J56_ -1/2 
2,295" 



)' 



C7 = - 



7,293' 
251,468^-1/2 



1/2 



\ 



(A24) 



C7 = 1(1 



C7 = - 



149,175' 
U\ ^ 34im^-l/2^ 

15,673 -1/2 

"7 <;CA /^ 1 



{\ + i2ZZ;,-l/2\ 

l,^ + 6,630" j' 



41 
'120 



(£5S-(fc)-i-lf.-'/^), 






_ 149 51 
~ 225 320 

8,113 9 

15,300 320 



346,687 
298,350 



41 
240' 



37 



5,850 ^ 80^'' 
15,362 , 11 



149,175 



480 



^, 



(A25) 
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with 



:^"', 


= c- 


~ 8,085 \^ " J^ 


^"', 


= C'!^ 


_ 1,024 1/2 
8,085" 


Ji"'^ 


= c- 


231 \^ 35" ; 



5. Mb = -1/2 



^rr ^ ^rr ^ ^(l_^-l/2)^ ^„ ^ ^„ ^ ''^^-''^ 

aiTT _ pTT _ 46-1/2 aiTS _ /^ts _ _4/'i _ 80t--l/2\ 

'^o ~ '-'o ~ 2i" ' "^-1 ~ '-'-I ~ zy ll'^ 7' 

^TT _ pTT -4/1 _ 69 -l/2~\ „„ _ r'-^s _ 148-1/2 

-^1 - <-l - 3l^ 28'^ j' A - <^0 - 231'^ ' 

^'-5^ -c^-i - -ili(i-"-'^^)' 

^SS _ f^SS _ 768 -1/2 

'-^-4 ~ ^^-4 ~ 2,695" 

<ziSS _ nSS _ 128/1 _ 73^-1 /2~\ 
^_3 - 0_3 - 77 ^^i -iqK j, 

rnSS _ /^SS _ 2,144 -1/2 rnvv _ ^w _ 992 

"^'-2 ~ '-'-2 ~ 2,695" ' -^-2 ~ ^ -2 ~ 8,085' 

J[S_l ^ c'_\ = -|(l-lf/.-l/2), ^- ^ ^vv ^ ^(l_ 82^-1/2)^ 

rnSS _ /^SS _ 1,002 -1/2 aiw _ /^w _ 788 ,,-1/2 

"^'O ~ '-'O ~ 385 " ' ^ ~ '-'O ~ 385" 

Cf =-li(£5S-'W-|f|"-'/^), C- = -|(55S-'«-^.-'/^), 

^SS _ pSS - l(] _ i864 -l/2'» raw _ ^vv _ 128-1/2 

-^3 - 1^3 - 4V 1,155'^ j' ^3 - C3 - -385/^ 

raTT _ prr _ .LOM/i _ -l/2^ 

■^-5 - '-'-5 - 8,085 \,^ " 7' 

^7T _ pTT _ _^12_ -1/2 

"^-4 ~ '-'-4 - 8,085" 

raTT _ prr _ .i^S/, _ 34 -l/2\ 

■^-3 ~ '-'-3 ~ 231 V 35" 7' 

^7T _ pIT _ _J68_ -1/2 

^i_2 - o_2 - 2,695" 

^- =c- = f(i-if.-'^^), 

<ZI7T _ ^TT _ 568 ,.-1/2 

^) ~ ^-'o ~ 105'^ 

^r - cr = -1(1 + 1.-'/^), 

J?[^^ = -f (55S"'('t) + TT - i""'^^) , 
C- = -f(55S-'W-i--'^^), 

^7T _ pTT _ i/l _ 512 -l/2^ 

^3 - O3 - 2 \,^ 385" 7' 



(A26) 



(A27) 



(A28) 



iRTT - 19 _ 3_ rt7T.s _ i„ _ 577 aySS _ 68,053 _ U 

-D - j7 4/1, -D _ ^y, ^g^, Z^ - 

8.278 3 „ (DTT _ 19 3 



- T^n, 



(DVV _ g,^/a i_ nrn — 12 _ J_ 

■" - 8,085 14^'' ^ ~ 1 r'- 



tlB 



(A29) 



The infra-red fields are dominated by a term proportional to k "^ ; in this case, this is A: . To ease the presentation, for this 
field we tabulate these coefficients separately from the others. 





J{% = C"3 


(JITT _ r>TT _ 1 

'^O - '-'O - 2' 


-^-2 = <^-2 


JiY = CY = -\\n{K), 


■^0 - '-'0 


rnTT _ /^TT _ 1 

^3 - O3 - 16, 


yi^ = CY 



flnW, 

3 

8' 

f, (A30) 



- C7 = -l(l-flnW), 

aiTS _ >^T5 1_ 

"^3 ~ *^3 ~ 32' 
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with 



and 



<:nS S _ /^S S 

-3 -3 16 V 3 '- 'J' aivv _ pvv _ _i n _ ]n(-J^,^^ 



-5 -5 8 ^ ^' <avv _ /^vv _ 1 ln/'l^^ 



<TiSS _ (^SS _ _i '•^ -2 ^-2 12' 

--^0 '^O 2' yjw - (^vv - _i 

^S5 _ ,-.ss _ 1 1 1 48 V 7 ^ ^y 



"^'3 ~ "-"S ~ 64 



yi!-; 


— 


C"', 


y[!'3 


= 


C'_\ 


^-2 


= 


C'^2 




= 




y[7 


= 


CY 



rnTT _ y^TT _ _! 

-^_4 - <-'_4 - 4^ 

j?[^3 = C"3 = |(l+41nW), 

^["2 = C"2 = 1^, (A32) 

"^0 ~ *^0 ~ 2' 



^- . Cr = -i(l-||lnW), 



_i_ 

3 - v^3 - 32> 



J?[o — C 



(A33) 



(A34) 



y[':'i = 3 (Li2(/t) + ln(/t) InW + I \n{K) + ^tt^) , 
C:\ = I (Li2(fe) - U2{k) + In(fe) InW + \ ln(^)) , 

y["j = I (Li2(/t) + ln(/t) ln(/f) - f InW + \ + j^jt^^ , 
C"i = I (Li2(fe) - Li2(^) + In(^) In(^) - | In(^) + i) , 

y[^^ = f (Li2(fe) + In(fe) In(^) + ^ In(^) + ^^n^ - ^) , 
C^f = I (Li2(fe) - Li2(^) + In(^) In(^) + ^^ In(^) - ^) , 

J?["; = 3 {hhik) + ln(/t) \n{K) + I ln(/^) + ^n^) , 
C"; = I (Li2(fe) - Li2(^) + In(fe) In(^) + \ ln(^)) , 

jr, - 9 (Li2(fe) + In(fe) In(^) - ^ In(^) +^^ + j^tt^) , 
C"^, = \ [U2{k) - \a2{k) + ln(fe)lnW - ^ InW + ^) . 

Appendix B: Analytic Power Spectra on Large Scales 

Only the largest scales are required for work with the CMB. On such scales the analytic solutions are significantly simplified; 
generating their Laurent series and keeping up to quadratic order in k, the analytic solutions for all hb e [-5/2, 3] reduce to the 
following. 

1. Wfl = 3 



V,Ak) « TiA^kl [\-k+ §k^) , Prsik) « nABkl {-\k+ ^^k^) , 

r^k) ^ nAskllf^ -k+ %e) , P„{k) « nAgkl {§-, - Ik + |f fe^) , (Bl) 
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where 



2. Mb = 5/2 



r,M) ^ nAek] (i^ - fe + iife2) ^ ^^^(^) ^ ^^^^3 (^ - 5^ + 1^2^ ^ (B2) 



riB 



n^W « nAekfa -k+ f-^e) , P„{k) ^ jiAekl (^ - f^ + f-,k^) , (B3) 

P„(k)^nABk^,{f,-lk+f^k^). 



4. A damped causal field 






c2f '^ 384 1 fV2e2f 

1 /l5V2;rerf(y2i)<;^-60f"2-8( 

~Ae^k+ ggo I f5/2g2f 



„ _ . ,3(_J_, 1 / 15 V2;r erf( y2i)g^-60g'^"-80f3^"+448f'^^-256f^/^ \ ,2 



_ uiilr _ J_t 1 / 1.365 V2^ erf( -^)e^-5,460f ■^^-7,285^^^^-12,992^^''^ +7,424f^'- ^ ,2 

/^ss — wA^fc^ I 5 -' ~ 7^'^ ~ 13,440 I ^5/2g2f I f!^ 

-P - 4 t3/'l4T - J^h 1 / 525 V2^ erf( V2f)e^-2,100f '/^-2,800g3/2-5,824f^/^+3,328f^/^ ^ ^^ 

"vv - ^AB/C^ljg7o 6g2f'f 6,720 I f5/2g2f I /^ 

_ Kiil&T 7 t 1 /735V2^erf(A/2i)e^-2,940f"^-3,920f^'2-10,304f"2+5,888f"2 



(B4) 



1 15 V2^erf(V2|)^Vf- 60^^/2 _ 30^9/2 _g4^i 1/2 
^0 = M ^WVf • (B5) 

5. Hb = 'ill 



P„(k) « nABkl {j-k+ l^k"-) , P,Ak) * nABkl {-\k + '^k^), 

Pss(k) « nABk', (t5 - ^ + ii^') ' 'Pvv(k) ^ ^A^fe^ (| _ 5^ + i«^2^ ^ (gg) 

^^(fe)«7rABfe3(56_7^^353^2~)_ 
6. Mb = 1 



^,,(^) ^ TtAb^^ (I _ fe + i|fe2~) ^ y,^^(^^ ^ ^^^^3 (I _ 5fc + |fe2~) ^ (B7) 

P„(k)^7TABkl{'^-lk+'^^k^) 
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7. Mb = 1/2 



Vrrik) = 


» ;rABfc3(i_fc_^fe2~)^ y,^^(^) ^ 7rAB^^(-^ 


Pssik) . 


. ;rABfc3(7_^_^^2~)^ ^^^(^^ ^ 7rABfc,3(|| 




nr(fe)«7rABfe3(2|_7^+ llfc2~) 



fe+3^fc'), (B8) 



8. Mb = 



r,Ak) « TrAfife^ (ff - fe + ^fc2) , P„{k) ^ nABkl\% -Ik- ^^e) , (B9) 

P„(k)^nABkl{'-§-lk+^^k^) 



Mb = -1/2 



Prrik) « TtAb^^ (2 - fe + I (in(ife) - iTT + fl) fe^) , 

!P„(fe) « nAsk', {-\k - (in(ife) - i;r + |) k^) , 

n,(fe) « ttAb^.^ (f - fe + H (in(ifc) - i^ + If ) fc') , (BIO) 

n.(fe) « ;rABfc3(28_5^+3(in(ife)-i;r+|f)fe2), 

^^(^) « ;rABfe,3(ff-|fe+f(ln(|)-i7r+f|)fe2). 



10. «B = -1 



!P^^(/t) » 7rABfe;?(4-5/t+|/t2)^ !P„(fe) ~ TrAs/t? (|/t - |f /t^) , 
^^(^)^7rABfe3(m_5^__i|^2~)^ 

11. hb = -3/2 

There is a logarithmic divergence for this field at low k, as part of the transition between the white noise nature of the stresses 
for hb > -3/2 and the scaling behaviour for hb < -3/2. We retain this logarithmic divergence and expand the remaining terms 
up to quadratic order in k; the dominant contribution as A: — > comes from the logarithmic term. 

Prrik) « nABk^ (-4 \n(\k)+f-27T-k+^k^), 

PrAk) « nABkllf-lk-'^k"-), 

P.dk) « TvABkl (-f In(ife) + If - M;r - fe + ^fe2) , (B12) 

P„ik) « nABk^ ("i ln(?^) + l " i^ " i^ " 3lj^') ' 
rAk) - nABkl (-i^ In(ife) + ^ - ff ;r - |fe - jife^) . 

On large scales, the ratio between the scalar cross-correlation and trace auto-correlation is 

Pr. 5 



Prr 91n(yt/4) 



0. 



The divergences in the other ratios cancel, and even though the white noise behaviour is broken, this transitionary field remains 

"causal". 
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12. ;ifl 



Pssik) ^ nAek] (f/r^^-i - f - ^) , V^k) ^ nABkl\ln^k-' - ff - f^) , (B13) 

VAk) ^ nAekl [In'k-' - '-^ - Ik) . 



The ratios between the spectra on the largest scales are therefore 

^ ' r5 _ ' 5.S' _ ^ ' VI' _ ^ 'TT _ ^ 
r TT f^TT ^ 'TT ^ 'TT ^ 

13. ub = -512 



r^k) ^ nABkl[^k-^-2-k), pTsik) * nAek^^k-^ - ^^k), 

Pss(k) ^ TiAekl (^r2 - ^ - fe) , VAk) * nAek] {fk'^ - i - g^) , (B14) 

'PAk)^7:Askl{Wk-'-^-lk). 



The ratios between the spectra on the largest scales are therefore 



23 n. 


Pss 


69 V,, 


69 r.-r 


23 


ATtt 


Vtt 


SOTtt 


188 ^T 


~ 17 



Appendix C: Discrepancy with Yamazaki et. al. 

In Yamazaki et. al. 1551 it was claimed [page 4] that "In almost all previous work the sum of the terms in brackets which 
include C in the k integral of Eq. (13) have been set to unity" and, in particular [page 7], that "These results are almost the 
same as in Brown and Crittenden [46] [61]. The constant ratios of all modes for ng > -1.5, however, are not unity because the 
terms including cosine factors [e.g. sum of terms within the bracket of Eq. (15)], is not unity as was assumed in the previous 
approximation." 

This is a misconception which we would like here to correct; in neither BC05, B06 nor the present work have we neglected 
any cosine factors. We briefly demonstrate this for the scalar trace auto-correlation (their equation (13)). The scalar trace auto- 
correlation, equation Jl2|l with Tj^ - \\^ '^ 1^} ^^^ tie rewritten as an integral across y, which is equal to Yamazaki et. al.'s C. 
Doing so, and assuming a power law spectrum, yields 

(r^) oc r ' fc'««+2 r' |k _ k'l"" (l + 1.?) dydk'. (CI) 

If we expand out // then this can be explicitly rewritten as 

(t2)oc I ;t'"«+2 r \\i-\i:\"''^U\+y^)k^-Akk'y + 2k'^)dyAk'. (C2) 

This is exactly Yamazaki et. al.'s equation (13); neither we nor they have neglected any terms. The "constant ratios" to which 
they refer are presumably the causal ratios presented in the previous Appendix, which are valid to the zeroth order in a Laurent 
expansion around k - 0. They are valid for all rig > -3/2 on sufficiently large scales; naturally, including the linear order in k 
will disrupt this. 

However, their statement (in their footnote [61]) that we did not fully consider the case with ng < -3/2 is, of course, entirely 
valid, as we merely considered these with realisations, an approach that is severely limited by an infra-red cut-off and by grid 
resolution, and did not attempt to integrate the expressions numerically, let alone analytically. 
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FIG. 1 : Sample realisations of the magnetic field and scalar pressures for a spectrum ng = -5/2 and side-length Zji,,, = 256. The magnetic field 
is on the back-right wall, the isotropic pressure on the back-left wall and the anisotropic pressure on the floor. 
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FIG. 2: Left: Probability distribution functions for the isotropic and anisotropic pressures. The j-axis is in units of the RMS amplitude of the 
field. Right: Skewnesses (black) and kurtoses (green) of the scalar pressures for iig e (-3, 2]. In the interests of clarity, negligible error bars 
have been suppressed on both plots. 
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FIG. 3: The stress spectra for a white noise magnetic field with iig = (left) and an inflationary field with iib = -5/2 (right). The infra-red 
damping of the realisations of the inflationary field atk < k^ is very apparent. 
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FIG. 4: 'Ptt (left) and P^s (right) for integer and half-integer n^ between ng = 3 (bottom, in black) and rig = -5/2 (top, in black). Also 
highlighted in black are rig = and ng = -3/2. 
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FIG. 5: Stress spectra for a causal magnetic field exponentially damped sXk -^ k^ (left), and an inflationary field with rig = -5/2 damped for 
k —^ (right). In the inflationary case numerical integrations are overlaid on the results for an undamped ng = -5/2 field. 
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FIG. 6: Left: Reconstructed power spectra for (top to bottom) iig = -2.9, «s = -2.75, n^ = -5/2, iig = 0, iig = 2 and the damped causal 
magnetic fields. Shaded regions denote the 1-cr errors. Right: Large-Scale CMB angular power spectra from magnetised tensor perturbations, 
normalised to /(/ + 1)C;|,^2 = !> for (top to bottom) iig > -3/2, rig = -3/2, rig = -2, rig = -5/2, ng = -2.75, ng = -2.9. Signals from the 
damped causal and IR-controlled inflationary fields are indistinguishable on these scales from those for ng = 2 and ng = -5/2 respectively. 



